Parameter Scaling in the Decoherent Quantum-Classical Transition for chaotic 

rf-SQUIDs 
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We numerically investigated the quantum-classical transition in rf-SQUID systems coupled to 
a dissipative environment. It is found that chaos emerges and the degree of chaos, the maximal 
Lyapunov exponent A m , exhibits non-monotonic behavior as a function of the coupling strength D. 
By measuring the proximity of quantum and classical evolution with the uncertainty of dynamics, 
we show that the uncertainty is a monotonic function of A m / D. In addition, the scaling holds in 
SQUID systems to a relatively smaller 7j e //, suggesting the universality for this scaling. 
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I. INTRODUCTION 

How classical behavior arises in a quantum mechanical 
system is one of the essential questions in quantum the- 
ory, and has long attracted intense interest. The quan- 
tum to classical transition (QCT), which has been well 
understood to be mainly induced by decoherence caused 
by the coupling with the environment,^ attains some 
progresses in recent years. It is proposed that the QCT 
is controlled by relevant parameters including the effec- 
tive Planck constant U e ff (i-e., the relative size of the 
Planck constant) , a measure of the coupling with the en- 
vironment D, and the Lyapunov exponent A, for chaotic 
systems.- By computing measures which directly reflect 
the "distance" between quantum and classical evolutions, 
it is shown that the distance is controlled by a composite 
parameter of the form ( = h a \ l3 D' y . Many efforts on 
investigating the coefficients a, j3, 7 have been made^£ 
in different systems such as the kicked harmonic oscilla- 
tor and the Duffing oscillator. However, in the previous 
systems, A is generally a constant. Therefore, the direct 
illustration of the effect of the Lyapunov exponent A on 
the computed distance is still open. 

In this article we try to explore the parameter scal- 
ing in QCT by using the system of the superconducting 
quantum interference device (SQUID). Rf-SQUID system 
has been demonstrated as a well controllable decoherent 
quantum system. Macroscopic quantum phenomena such 
as resonant tunneling^ and level quantization— and quan- 
tum superposition^ have been reported. On the other 
hand, the strong coupling between the SQUID and the 
environment can introduce chaos. As early as 1983, the 
chaotic behavior of the SQUID treated as a semi-classical 
model had been founds Recently, a research shows that a 
three-junction SQUID can be used to study the dynamics 
of quantum chaosJ^ Such works motivate us to study the 
chaotic behavior of SQUID under decoherence induced by 
environment, which enables us to directly demonstrate 
the effect of the Lyapunov exponent on QCT. 

This article is organized as follows. In Sec. II we numer- 
ically investigate the chaotic dynamics of SQUID with 
coupling to an external environment, and it is shown that 
the maximal Lyapunov exponent A m , which quantifies 



the chaotic degree of SQUID, is non-monotonic as a func- 
tion of D, a measure of the coupling. Thus we can say in 
some regimes of D, the chaos of SQUID is suppressed by 
the decoherence induced by environment^. In Sec. Ill we 
use the uncertainty of dynamics as the distance between 
quantum and classical evolutions, and show that the un- 
certainty behaves rightly, even in the chaos suppressed 
region, as a monotonic function of \ m /D. To the best of 
our knowledge, this is the first direct demonstration of 
the scaling relation since it was proposed 3 . 



II. CHAOTIC DYNAMICS OF SQUID 

The rf-SQUID system considered here consists of 
a large superconducting loop interrupted by a single 
Josephson junction with a critical current I c . Under 
the driving of a external flux <j) ex {t) with the form of 
4>ex(0) cos(u!dt) (where ex (O) and ujd respectively denote 
the driving amplitude and driving frequency), the Hamil- 
tonian for the SQUID system can be given as 



Hd ~2C 



(0-<Mt)) 2 

2L 



2tt 



■ cos(2n<j>/<j> ), (1) 



where C is the junction capacitance, L is the rf-SQUID 
inductance and </>o = h/2e denotes the superconduct- 
ing flux quantum. The magnetic flux threading the rf- 
SQUID cf> and the total charge on the capacitor q are 
the conjugate variables of the system with the imposed 
commutation relation [</>, q] = ih. 

We can rewrite this Hamiltonian into a dimensionless 



Q 2 , ($-<M*)) 2 



■ cos( 



2e 



2u) e y/haj C 



*), (2) 



Hd = t 

in which u>q 

Q = y/l/huJoCq, $ = \Juj$C/h(j) satisfy the commu- 
tation relation [$, Q] = i. 

Since no chaos can be seen in the dynamics of iso- 
lated quantum systems^ to study the chaotic behaviors 
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(a): D=0.25 " 
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(b): D=0.35 
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FIG. 1: Poincare sections for D = 0.25, 0.35, 0.45, from top to 
bottom. From middle panel we can see that points are largely 
confined in three regions, which indicates a non-monotonic 
transition of chaos. 



of the SQUID system, we couple the system to a dis- 
sipated environment in the Markovian limit. We adopt 
the quantum state diffusion (QSD) 14 approach which is 



15.16.17 



to 



widely used in studying open quantum systems 
describe the evolution of this coupled system. The QSD 
equation for the evolution of the state vector reads 



3^3 



(L]){Lj)) \i>)dt + - (4-))|^)^(3) 



where H is the system Hamiltonian and Lj are the Lind- 
blad operators representing the coupling with the envi- 
ronment. d£j are independent complex differential Gaus- 
sian random variables satisfying M(d£j) = M(d^id£,j) = 
0, M(d£*d£j) = Sijdt (where M denotes the ensem- 
ble mean). For the SQUID system considered here, we 
have H and L for Equation ((3]) as H = Hd + Hr, 
L = \/15(Q + iQ), where Hjj is shown in Equation |(3J), 
Hr = — ($Q + Q<f>)i£ii& is a damping term added to re- 
cover the correct equation of motion in the classical limit, 
and D is the strength of the coupling with the environ- 
ment mentioned in the beginning. 

Using the powerful QSD library^ we numerically 



solve the Equation ([3]) and investigate the change in 
the dynamics of the SQUID system when increasing the 
strength of dissipation. A typical set of SQUID pa- 
rameters is selected here, C = Q.lpF, L = 'SQQpH, 
I c = 2.2,uA, ui d = 1.14w , <Pex(0) = 0.268400, which in- 
sures the action of this system is small enough compared 
with fixed hJ£ Then we examine 28 different values of D 
from slightly dissipated (D = 0.23) to heavily dissipated 
regime (D — 1) in our calculation, during which we have 
the same initial state \ip(t = 0)) = |%/2(($) + i(Q)) = 
(0.877 — 0.566z))-the coherent state-and same realiza- 
tion of generating the random numbers. The quantum 
Poincare sections, which each comprises of 500 points 
taken at a fixed phase of the external driving once a driv- 
ing period, are shown for three representative values of 
D in Fig.l(a)-l(c). It can be clearly seen in Fig. 1(a) that 
points forms a uniformly stretched Poincare profile in the 
phase space which indicates "chaos" for D — 0.25. How- 
ever, for D — 0.35 most of points are confined in three 
relatively small regions as shown in Fig. 1(b), which indi- 
cates the suppression of chaos. Then the Poincare pro- 
file similar to the one in Fig. 1(a) is recovered in Fig. 1(c) 
when D is increased to 0.45. Some non-monotonic anal- 
ogous phenomena have been studied in classical chaotic 



systems 



11.19 



and a qualitative explanation has been pro- 



posed there. If the chaotic attractors are narrowly and 
non-uniformly distributed in phase space, the fluctuation 
induced by dissipation may cause the neighboring tra- 
jectory jump over it, which results in the suppression of 
chaos. While further increasing dissipation intensity, the 
structure of the chaotic attractor may be modified and 
thus spread wider than before. Therefore the system be- 
comes chaotic again. Since (($(t)), (Q(t)}) form classical- 
like trajectories in our calculation, we expect that the 
explanation is also valid for the suppression of chaos in 
quantum region. 
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FIG. 2: Maximal Lyapunov exponent A m versus D.The dis- 
tinctive dip rightly attests the occurrence of suppression of 
chaos. 



To describe this transition of chaos quantitatively, we 



3 










0.95 




(b) 






0.90 










- 


■ 




0.85 












0.9 ~ 


- 


■ 




0.80 












0.8- 


- 






0.75 


<,« 0.7- 














■ i 

■ i 




0.70 


0.6- 










0.5- 








0.65 


0.7 








0.60 




■ ■ ■ 






0.55 




-3.0 -2.5 


-2.0 -1.5 -1.0 


-0.5 





Iog 10 (i/D) 



iog 10 (VD) 



FIG. 3: Averaged uncertainty A a as a function of (a) D and 
(b) a composite parameter \ m /D. The monotonic increase 
of A a as a function of \ m /D in (b) demonstrated the scaling 
law. 



calculate the maximal Lyapunov exponent X m for a time 
series-the expectation value of the magnetic flux 
at each value of D. The calculation is based on the 
method and program o 20 i 21 which are specifically designed 
for the analysis of nonlinear time series. With carefully 
chosen parameters as the delay time d — 3, the embed- 
ding dimension m = 3 and the scaling length s = 1.4% for 
the calculation to best meet the requirements in Ref.18, 
the sufficient convergency of the Lyapunov exponent is 
guaranteed. The result is shown in Fig. 2, in which the 
graph of X m versus D has a distinctive dip in a approxi- 
mate region of D = 0.25 ~ 0.45, indicating the suppres- 
sion of chaos. We also repeat the whole calculation above 
in some different realization of random numbers with the 
SQUID parameters and the initial state fixed, and find 
the curves are quite analogous to the one in Fig. 2. 



III. EFFECT OF MAXIMAL LYAPUNOV 
EXPONENT ON QCT 

With the non-monotonic relationship between maxi- 
mal Lyapunov exponent X m and the strength of the cou- 
pling with the environment D, we can directly inves- 
tigate the effect of X m on QCT. To measure the "dis- 
tance "between quantum and classical evolution, we use 
the well known quantity-the uncertainty of dynamics 

A : 



'(($ - (&)) 2 )y/({Q - (<3)) 2 ), which is simple for 
calculation and adequate to describe the QCT. Accord- 
ing to the commutation relation Q] = i, it follows 
that A > 0.5. By solving Equation ([3]) with same cal- 
culating parameters as in Sec. II, we get a time series of 
the uncertainty A(t) at each value of D. After averaging 
each series of A(t) over a reasonably long time (> 100 
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FIG. 4: (Color online) Averaged uncertainty A a as a function 
of D and h e ff. The parameters for the system with largest 
effective Planck constant h e ff = 2.6 has been shown in text. 
Other h e f f and corresponding sets of parameters are listed in 
Tabled. 



TABLE I: heff and corresponding parameters 





I D <jxA) L( P H) C(pF) 


UJ d (uj ) <j) 


^(0)(^ ) 


1 


4.6 100 3.95 


0.65 


0.081 


1.2 


3.35 150 2.16 


0.71 


0.1041 


1.4 


2.67 200 1.29 


0.78 


0.1273 


1.9 


2.46 250 0.36 


0.99 


0.1851 



periods of the external driving), we obtained the curve 
of the averaged uncertainty A a versus D and showed in 
Fig. 3 (a), where D has the same sequence of values as in 
Fig. 2. It can be clearly seen that in Fig. 3(a) a obvious 
dip emerges in the very regime where chaos is suppressed 
by the dissipation, which implies QCT directly depends 
on the degree of chaos. Motivated by this, we combine 
X m and D with the form of X m /D which is inferred in 
Ref.2 and look into the relationship between A a and such 
composite single parameter. Shown Fig. 3(b) is an exam- 
ple of A a vs. X m /D. One can find that the dip is rubbed 
out and A a approximately shows a monotonic increasing 
in X m /D with two distinct regimes of small and large 
increasing rates £ Therefore we demonstrate the scaling 
between A m and D holds over a considerable range in 
A a . It is noticed that the points which lie in the dip 
in Fig. 3(a) spread slightly around the curve in Fig. 3(b). 
We conjecture this spread could be mainly attributed to 
the calculating error— of X m which is induced by the 
inevitable quantum noise added into the trajectory of 
(($(t)), (Q(t))), especially when chaos is suppressed and 
the value of X m is comparatively small. 

Now we examine this scaling law for the SQUID system 
with a smaller effective Planck constant h e ff. To obtain 
a smaller ?i e //, it is not straightforward for the SQUID 
system to directly manipulate the value of Ti ^Instead, 
we enlarge the action of the SQUID system simply by 
changing parameters in the Hamiltonian; the larger the 
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comparison. Fig. 4 shows the averaged uncertainty A a as 
a function of D, h e ff. For each 7i e //,a distinct dip ex- 
ists as expected in the region where chaos is suppressed 
by the dissipation of environment. Fig. 5 shows the same 
data plotted as a function of \ m j D, in which, the behav- 
ior of A Q for each h e f / is considerably the same, which 
demonstrates the scaling between A m and D is still valid 
for a system with relatively small % e ff- For clarity, we 
separately show the curve with fr e ff = 1 hi the inset of 
Fig. 5. Since a larger action is helpful to undermine the 
effect of quantum noise, more accurate A m can be gained 
for the system with smaller h e f f , which, is reflected in the 
lack of noticeable spread around the curve in the inset. 

We also chose some different random numbers gener- 
ator to repeat the calculation for SQUID systems with 
different h e ff> and succeed in getting same qualitative 
conclusions as discussed above. 



FIG. 5: (Color online) A a versus a composite parameter 
Xm/D for different effective Plank constant. It is shown that 
the scaling law holds for systems with different h e ff. Inset: 
The curve with Ti e f / = 1 is shown separately. 

action the smaller K e f / , and vice versa^ By deliberately 
selecting the set of parameters including I c , L, C, u>d 
and 4> ex (0), we can enlarge the action and maintain the 
chaotic dynamics of the system at the same time. The 
values of these parameters are not difficult to modulate 
for a realistic SQUID system where I c could become con- 
trollable by replacing the single Josephson junction with 
a small loop (dc SQUID) which contains two identical 
Josephson junctions, - C and L are both under the upper 
realistic limit of typical Josephson junctions. We select 
four sets of parameters for the SQUID systems each of 
which has a smaller h e f t compared with the foregoing 
system's. Assuming the smallest % e ff is equal to 1 and 
comparing the actions of the systems which are measured 
with the system size^ we approximately gain the value 
of other effective Planck constants as follow, 1.2, 1.4, 1.9, 
2.6, where 2.6 is the value of the foregoing system's h e ff. 
Then we apply the same calculating procedures to these 
systems, and the results are shown in Fig. 4 and Fig. 5 
which also include the data of the foregoing system for 



IV. CONCLUSION 

In summary, we investigated QCT in chaotic rf- 
SQUIDs. The suppression of chaos induced by environ- 
ment dissipation was observed in quantum regime. It 
is found that the quantum to classical transition in the 
presence of a dissipated environment is governed by a 
composite parameter X m /D. It could be expected the 
scaling law between A m and D would holds over a wide 
range of h e ff. However, to generalize this scaling to the 
one involving h e ff, A m and D and to reveal the coeffi- 
cients between them are still open questions needed to 
explore. 
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